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ANALYTICAL AND EXPERIMENTAL INVESTIGATION OF TEMPERATURE RECOVERY 
FACTORS FOR FULLY DEVELOPED FLOW OF AIR IN A TUBE 
By R. G-. Deissler, W- F. Welland, and W. E. Lowdermilk 


SUMMARY 

An analysis was made for predicting temperature recovery factors for 
fully developed flow in a tube. Most of the attention was confined to 
turbulent flow. Some qualitative results were obtained for laminar flow 
by setting the eddy diffusivity in the equation for turbulent flow equal 
to zero and using the incompressible parabolic velocity profile for lam- 
inar flow. For zero Mach number the laminar flow results were exact. 
Radial variation of properties was neglected in most of the calculations. 
The effect of wall temperature gradient along the tube was negligible for 
turbulent flow below Mach numbers of 0.9 and 0.98 for Reynolds numbers of 
20,000 and 390,000, respectively. For laminar flow the effect became im- 
portant at much lower Mach numbers. 

Recovery factors were obtained experimentally for a range of Reyn- 
olds. number from 630 to 30,000. Additional previously unpublished data 
are presented for Reynolds numbers up to 650,000. The results indicate 
that in the turbulent flow region the recovery factor is approximately 
independent of Reynolds number. In the transition region for Reynolds 
numbers between 2000 and 3000 the recovery factor is reduced abruptly to 
a value lower than that obtained for the turbulent flow region. 

Comparison of the predicted recovery factor with experimental data 
indicated that the effective value of ratio of eddy diffusivities for 
heat tr ans fer to momentum transfer varies from 1 at a Reynolds number 
of 5000 to 1.09 at a Reynolds number of 400,000. The ratio of eddy 
diffusivities was also obtained from measurements of heat-transfer coef- 
ficients and an analysis . Comparison of predicted with measured heat- 
transfer coefficients indicates effective eddy diffusivity ratios slightly 
less than those calculated from recovery factors. 

The analysis substantiated the experimental result that correlations 
for heat transfer without frictional heating can be applied to the case 
with frictional heating hy basing the heat- transfer coefficient on the 
difference between the wall temperature and the adiabatic wall tenqperature 
rather than on the difference between the wall and the fluid hulk 
temperature. 
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INTRODUCTION 

Most of the work on temperature recovery factors has been confined 
to supersonic flow over external surf aces (ref. l) . Frictional heating 
at high supersonic Mach numbers, together with material limitations, nec- 
essitates accurate calculations of surface temperatures in this region. 

A knowledge of recovery factors is also necessary for calculating 
convective heat transfer for flow in tubes at high subsonic Mach numbers 
and small temperature differences. This case was studied experimentally 
and reported in reference 2. The work was extended to supersonic flow 
in the study of reference 3. It is shown in reference 2 and by the anal- 
ysis herein that, for turbulent flow, low-velocity heat-transfer correla- 
tions can be used for higher velocities if the heat- transfer coefficient 
is based on the difference between the wall and the adiabatic wall tem- 
peratures rather than on the difference between the wall and the fluid 
bulk temperatures. Calculation of the adiabatic wall temperature, or the 
temperature the wall would assume with no heat transfer, depends on a 
knowledge of the recovery factor. 

The present experimental work covers a range of Reynolds numbers 
from 640 to 30,000. The values reported in reference 2 were obtained 
for a Reynolds number of approximately 30,000. Previously unpublished 
values obtained in connection with the work in reference 4 were measured 
for Reynolds numbers from 250,000 to 650,000 and are reported herein. 

Very little analytical work has been done on fully developed recov- 
ery factors, although some analysis appears in reference 5. (Fully de- 
veloped is herein taken to mean that the recovery factor is independent 
of distance from the tube inlet and that the velocity and total- 
temperature profiles are similar at various distances from the tube in- 
let.) In the present report, most of the attention is confined to fully 
developed subsonic turbulent flow. Some qualitative results are given 
for laminar flow by dropping the eddy diffusivity from the equations for 
turbulent flow. 

In the turbulent case the recovery factor is quite sensitive to the 
ratio of eddy diffusivities for heat transfer to momentum transfer. The 
mean effective eddy diffusivity ratio at a crass section of the tube can 
be estimated by comparing the present analysis with the measured values 
of recovery factors reported herein. 


SYMBOLS 


C constant of integration 

c specific heat of fluid at constant pressure 
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D 

h 

k 

n 

P 

q. 

% 

R 

r 

r 0 

T 



^0,a 

t 

*a 

% 

^b,a 


u 

^b 


Inside tube diameter 

heat-transfer coefficient for no frictional heating, q^/(tQ - t^) 
heat-transfer coefficient with frictional heating, q^/ (tg - tQ a ) 
thermal conductivity of fluid 
constant, 0.124 
static pressure 

rate of heat transfer per un i t area toward tube center 

rate of heat transfer per unit area from, inside wall toward 
center 

gas constant 
radius 

inside tube radius 

2 

total temperature, t + - 5 —-, deg abs 

^p 

bulk or average total temperature of fluid at cross section of 
tube, deg abs 

bulk or average total temperature of fluid at cross section of 
tube with no heat transfer, deg abs 

wall temperature , deg abs 

adiabatic wall temperature, deg abs 

static temperature , deg abs 

static temperature with no heat transfer, deg abs 

b ulk or average static temperature of fluid at cross section of 
tube, deg abs 

bulk or average static temperature of fluid at cross section of 
tube with no heat transfer, deg abs 

velocity para ll el to axis of tube 

bulk or average velocity at cross section of tube 
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velocity in radial direction 


V 

velocity in negative radial direction 

X 

distance along tube 


y 

distance from tube wall 


r 

ratio of specific heats, taken 

as 1.40 for air 

6 

eddy diffusivity for momentum 


6 h 

eddy diffusivity for heat 


X 

Karman constant, 0.36 


P 

viscosity 


V 

kinematic viscosity, p/p 


p 

density 



bulk or average density 


T 

shear stress in fluid 


T o 

shear stress at wall 


Dimensionless Parameters: 



eddy diffusivity ratio, e^/ e 


f 

D(dp/dx) fr 

friction factor, 5 — — = 

2 P*b 

u b 

Mach number. 

. fo 

M 


Nu 

Nusselt number for no frictional heating, hD/k 

Nu. 

0,a 

Nusselt number with frictional heating, h^gD/k 

Pr 

Prandtl number, c^p/k 


Re 

Reynolds number, pu^D/p 



to* 


T62t. 
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T+' 

< 

*b 

u + 


+ 

u b 


U.n 


4 * 


<p 


tube radius parameter. 


W pr o 


total- temperature parameter for adiabatic case. 


2(T n - T)c p 


total bulk temperature parameter for adiabatic case, 
2(T q - T b )c.pP 


0 


bulk temperature parameter. 


(*0 - t b) c p T 0 

t o/ P 


adiabatic wall temperature parameter. 


fro - Wyo 

9oV T o/ p 


velocity parameter, u/ -/ T q/ p 
bulk velocity parameter, u^/i/tq/p 

value of u + at y^ 

wall distance parameter, (V T o/ P y)/v 

value of y + at intersection of curves for flow close to wall 
and flow at distance from wall 


recovery factor , 1 


*b - *0 
U b/ 2c p 


Subscript : 

fr on friction-pressure gradient 


AHALYSIS 

In order to calculate the recovery factors for flow through a tube, 
the temperature distribution in the tube will first be obtained from the 
energy equation. The energy equation for fully developed turbulent flow 
in a tube is derived in appendix A and can be written as 

( r o - y> c P pu H = b [ (r ° ' y) ( k iy " pc p v,t ’ + ^ tp • pu u ’ v ')] (1) 
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where T is the total temperature and t the static temperature. The 
bars denote time averages, and the primes indicate fluctuating compo- 
nents . Equation (l) is the same as equation (A7) in appendix A if the 
bans over the time-averaged velocities, temperatures , and properties are 
dropped. 


By introducing the eddy diffusivities for momentum and heat trans- 
fer, defined by 


u'v' 



t'v’ = 



and setting 

t = (n + pe) ~ 


( 2 ) 


equation (l) becomes 

( r 0 “ y> c pP U Is = ^ j( r 0 - y) [(* + pc p ^e) + u-rj | (3) 

where ^ is the ratio of eddy dlffusivities for heat transfer to momen- 
tum transfer. 


The present section concerns the case where the wall is insulated, 
and hence 


St 

& 


= 0 


y=o 


(4) 


If equation (3) is multiplied through by dy and integrated between 
the wall and the center of the tube, there results 


| f pu3? ( r o - ytoj = f d |( r o - y)[( k 


, \ St , 

+ pc p <3,e) + UT 


! 


( 5 ) 


where pu was regarded as independent of x from continuity. The 
right side of the equation is zero because the expression in brackets 
is zero at the insulated wall and at the center of the passage. The 
bulk temperature is defined 
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m = 3 
i b- — 


f ° pul(r 0 - y)dy 

•A 


/ ± 0 

pu(r Q - y)dy 


( 6 ) 


for constant specific heat. The numerator and deno min ator in equation 
(6) are independent of x "by equation (5) and continuity. Therefore 


ar b 

dx 


= 0 


(7) 


This result is used in the next paragraph. Equation (7) can, of course, 
he obtained by other methods, but it is of interest that it follows from 
equation (3). Although dT^/dx = 0 by equation (7), dT/dx f 0 at some 
radii, so that the left side of equation (l) will be retained and an es- 
timate made of its importance. 

It is assumed herein that the total- temperature profile is fully 
developed, or 


m m" = /(y) (S) 

0 x b 

Differentiation of equation (8) with respect to 

dT ^0 / ^0 “ T \ 

dx - dx \J- - T 0 - T b ) 

where equation (7) was used. 


Dimensionless Form of Energy Equation 

In the remainder of the analysis the effect of radial variation of 
properties is neglected. It is shown in reference 6 that the effect of 
radial . variation of properties on the turbulent velocity profile is neg- 
ligible for adiabatic flow at subsonic Mach n umb ers . It is reasonable 
to assume that the same result holds for turbulent temperature profiles 
and recovery factors. 

If equation (9) is substituted into equation (3) and dimensionless 
quantities are introduced, the following equation results 


x gives 

(9) 
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where du/ dy 



was eliminated "by equation (2) and 


A s 


STq 

~ST V 
T 0 VV^ 


( 11 ) 


In obtaining equation (10) , the static temperature was eliminated 
by using the relation 


T = t + 


u 

2c, 


( 12 ) 


An estimate of A is made in appendix; B by using one-dimensional 
flow theory. One -dimensional flow is a reasonable assumption for tur- 
bulent flow because the velocity and temperature profiles for that case 
tend to be uniform over most of the cross section. The resulting expres- 
sion for A is 


A 


2*b' r M 2 
rjuf 1 - M 2 


(13) 


For the quantity t/ Tq 
relation 


in equation (lO), 


the incompressible 



(14) 


will be used. It Is shown in reference 6 that the velocity profile in 
turbulent flow is insensitive to the way in which the shear stress varies 
with radius. It would be expected that the same result would hold for 
recovery factors. 


The use of equation (14) for compressible turbulent flow can also 
be justified by consideration of the momentum equation (eq. (Al) ) . Equa- 
tion (Al) can be used directly for fully developed turbulent flow by re- 
placing p by p + pe and letting u r = 0. For turbulent flow, where 

the velocity profile tends to be flat, u du/Bx « u^ Su^/cbc. If the 
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pressure gradient is separated into the momentum pressure gradient 
-pu^ cki^/cbc and the friction pressure gradient -ZXqJxq, and equation 

(2) is introduced into equation (Al ) , equation (14) results on 
integration. 

Substituting equations (13) and (14) into (10) and integrating 
twice gives 


+ 



Expressions for Eddy Diffusivity 

The eddy diffusivity e in equation (15) is obtained from the fol- 
lowing expressions, which are experimentally verified in references 4 and 
7 : For flow close to the wall (y+ < 26) 

e = n 2 uy(l - e- n *Wv) ( 18 ) 

For flow in the region away from the wall (y + > 26) the Karman relation 

e = x 2 (du/dy) 5 ( 19 ) 

(d 2 u/«ay 2 ) 2 

is used, where n and x are experimental constants having the values 
0.124 and 0.36, respectively. 
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By integration of equation (2) for the region away from the wall 
with viscous shear stress neglected, and using equations (14) and (19), 
von Karman obtained the following dimensionless equation (ref. 8) : 



Equation (19) becomes, on substitution of the first and second deriva- 
tives from equation (20), 



for flow in the region away from the wall. 

For flow close to the wall, equation (18) can be written in dimen- 
sionless form as 

| = nV/(l-«- A V) (22) 

where the values of u + are obtained from the generalized velocity dis- 
tribution in figure 1. 

Values of T+' as a function of y + for various values of r£ 
and Mach number can now be calculated from equation (15) . 


Temperature Recovery Factors 


The definition 
(ref. 2), and which 


of recovery factor usually used in experimental work 
is used here, can be written 


<? = 1 


^b - T 0 
u|/2c p 


An alternate definition which is sometimes used is 


(23) 


t 0 ~ ^ 

v£/ 2c v 


(24) 
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However, equations (23) and (24) are strictly equivalent only for a uni- 
form velocity distribution, inasmuch as 


®b = *b + 2^ (®) 

where a is a factor accounting for the nonuniformity of the velocity 
profile, having a value of 1 for a uniform profile. Inasmuch as total, 
rather than static, bulk temperatures are measured in the present experi- 
ments, the use of equation (23) is to be preferred. 


Equation (23) can be written in terms of 
ties as 


readily calculated quanti- 


(P = 1 



(26) 


where and u^, for a given r£ and M, are obtained from equations 

(15), (16), and (l7) and the u + against y + plot in figure 1, obtained 
from reference 7. If it is desired to obtain cp as a function of the 
usual Reynolds number rather than of r§, the Reynolds number can be cal- 
culated from 


Re = 2Ub r 0 


(27) 


External Heat Transfer 

For calculating the external heat transfer for flow thro ugh a tube 
at high velocities, the energy equation (3) can be written as 

( r o - li = " I? r o - y)i] (28) 

where 

q = -(k + pc <Ze) - ut (29) 

P dy 

For the case of no external heat transfer from the wall, it is shown in 
the section Effect of Hall Temperature Gradient on Recovery Factor that 
the left side of equation (3) can be neglected except at very high sub- 
sonic Mach numbers. Equation (28) then shows that q « 0 at aJ 1 radii. 
For no external heat transfer equation (29) therefore becomes 


0 = 


-Os + Pdp&e) 




UT 


(30) 
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where ut Is assumed to he unaffected by heat transfer (small heat 
transfer) . Subtracting equation (30) from equation (29) gives 


q. = 

Integrating equation (31) 
t a - t + t 


-(k + pc p <2.e) 


d'(t - t a ) 

fly 


from the wall to y gives 

•y 


0 “ t 0 


..-I 

Jo 


(k + pc p a-e) 


fly 


(31) 


(32) 


Multiplying equation (32) by (r Q 
the tube center, and dividing by 


- y)u dy, integrating from the wall to 



give 


%,&, " ^ + ^ “ ^0,a 


r r o 

f 7 1 


/ u(:r 0 - y) 

/ (k + pc O-e) ^ 

fly 

Jq 

L/o J 




- y)fly 


(33) 


Setting a , or T b = T b a> since 1:116 k eafc transfer is assumed to 

be small and thus the effect of heat transfer on the velocity profile can 
be neglected, results in 



^0,a 




+ pc p <2,e) 


fly 


fly 



- y)fly 


(34) 


The quantity t^ was set equal to a because the wall temperature is 
co mp ared to the adiabatic wall temperature for the same bulk temperature. 

In dimensionless form equation (34) becomes 
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q/qQ dy + 

1 , ae 
Pr + v . 



dy + 


(35) 


The exact variation of q/q q in equation (35) is not important for 

turbulent flow (ref. 4), except possibly at very low Reynolds numbers. 

For 3 t/Sx independent of y, equation (28) can be integrated first, 
from the wall to a point in the fluid, and then from the wall to the cen- 
ter of the tube. Eliminating cfC/chc between the two equations and con- 
verting the result to dimensionless form give 


= 

“30 


-+ 



r v+ 

t $ *4 u+(r o ■ y+)ay+ 

r o ■ y+ r r o 

J 0 u+ ( r Q - y + )dy + 


(36) 


Equation (36) is the same as the equation for no frictional heating ob- 
tained in reference 8. Note that q as used herein contains a contribu- 
tion from dissipation (eq. (29)). 


The Nusselt number is given by 


Nu 0,a 


2rJPr 

^0,a 


(37) 


Equation (37) follows from the definitions of the various quantities 
involved. 


If in equation (29) the term ut had been neglected, that is, if 
frictional heating had been neglected in the analysis, t* in equation 

U y EL 

(35) would have been replaced by t^, and the equation for Nusselt number 
would become 


Nu 




(38) 


Thus, the result for no frictional heating is identical to that for fric- 
tional heating if the bulk temperature in the definition of the heat- 
transfer coefficient is replaced by the adiabatic wall temperature . This 
is the result which was found experimentally in reference 2. 
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APPARATUS AND PROCEDURE 

The present apparatus was designed to .determine the effect on re- 
covery factor of Reynolds numbers less than 30,000. Air at subatmos- 
pheric pressure was passed through a thin-walled tube test section. 


Arrangement of Apparatus 

The general arrangement of the apparatus is shown in figure 2. Air 
at room temperature and pressure first enters a temperature stabilizing 
box with a volume of 2 cubic feet packed with steel wool. Then it flows 
in turn through a throttling valve, a rotameter, an inlet ml xi ng tank in 
which the total temperature and pressure of the air are measured, the 
test section, an outlet mixing tank in which the outlet total temperature 
is measured, a second throttling valve, an exhaust tank, and the labora- 
tory altitude exhaust system. The density and velocity in the test sec- 
tion were controlled by the throttl ing valves. 

The mixing tanks were thermally insulated by a layer of glass wool 
enclosed in aluminum foil. The test section was thermally insulated by 
a layer of glass wool enclosed in a thin steel cylinder 3 inches in 
diameter. The entire apparatus was enclosed in a tent to minimize the 
effects of fluctuating room air temperature and currents produced by the 
room air conditioner. 


Test Section 

A diagram of the test section is shown in figure 2. The test sec- 
tion was the middle portion of an Inconel tube 51 inches long with an 
inside diameter of 0.222 inch and a wall thickness of 0.014 inch. The 
distance between the first and last pressure tap was 44.4 inches and com- 
prised the test section proper; thus the approach and exit sections were 
short (3.3 in.). The test-section tubing was fastened to the mixing 
tanks by means of 1-inch cubical Lucite blocks which minimized the con- 
duction of heat at each end of the test section. 

Static pressure was measured at the entrance and exit of the test 
section by means of pressure taps located in Lucite blocks. 

The thermocouples consisted of 32-gage .iron- const ant an wire, first 
butt-welded to form the junction and then spot-welded to the test section. 
Location of the thermocouples is shown in figure 2. The wall temperature 
thermocouples were connected differentially 'with the inlet total-air tem- 
perature thermocouples, and the differential temperature was measured by 
means of a self-balancing potentiometer modified to give full-scale de- 
flection at 1 millivolt. 
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Procedure 

Before each run the reading on the potentiometer for no flow through 
the test section was checked to obtain the zero reading. Then the air- 
flow rate and pressure level in the test section were set at the desired 
values. For high flow rates temperature equilibrium was usually reached 
in a few minutes , and all readings were recorded. For low flow rates 
temperature equilibrium was difficult to discern, and readings were re- 
corded shortly after the difference between the inlet total temperature 
and the wall temperature had reached a maximum value. Readings were 
thereafter recorded whenever the temperature difference reached a max- 
imum or minimum value for a 2- to 3-hour period. 

Measured values of static pressure at the inlet of the test section, 
static -pres sure drop along the test section, airflow rate, inlet total 
t emp erature of the air measured in the inlet mixing tank, and temperature 
difference between the inlet total temperature and the adiabatic wall 
temperature at the exit of the test section are given in table I. 


RESULTS AM) DISCUSSION 
Experimental Results 

Values of temperature recovery factor as defined by equation (23) , 
average friction factor, and Reynolds number obtained in the present in- 
vestigation are given in table I. The recovery factor is based on the 
measured temperature difference between the total air temperature in the 
inlet tank and the wall temperature at the exit of the test section. For 
airflow rates corresponding to Reynolds numbers above 3000, the measured 
temperature difference was affected slightly by fluctuations in the room 
air temperature of 1° or 2° F. The total air temperature thermocouple 
responded faster to changes in the air temperature than the wall temper- 
ature of the test section; hence, an increase in room air temperature re- 
sulted in a maximum difference between the total air temperature and the 
wall temperature, whereas a decrease in the room air temperature resulted 
in a nHnlmum temperature difference. The room air temperature variations 
were caused by the cyclic operation of the room air-conditioning unit 
and were usually repetitive at half-hour intervals. The time required 
to reach equilibrium was only a few minutes at the higher flow rates, and 
hence the minimum and maximum values of temperature difference deviated 
only slightly from the time-averaged value, and the recovery factor can 
be indicated by the numerical average of the minimum and maximum values 
(e.g. , the last six runs). 

A strictly adiabatic flow condition was difficult to obtain in the 
experimental work, since a temperature difference existed between the sur- 
rounding air in the room and the wall of the test section. Therefore, 
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calculations were made to determine the effect on recovery factor of 
heat transferred to the air in the test section from the surroundings, 
even though there was no indication "by measurements that the total air 
temperature increased between the inlet andexit of the test section. 
The calculations indicated that this effect on recovery factor was 
negligible . 

The estimated maximum possible error in calculating the recovery 
factor is as follows: 


Reynolds number 

Error, 

percent 

Above 16,000 

±2 

4,000 to 16,000 

±6 


Below a Reynolds number of 4,000 the maximum possible error ranges from 
±8 to ±26 percent, and 70 percent of these data have an error of less 
than ±20 percent. Approximately two-thirds_ of the total error can be 
attributed to the combined errors of measuring T^ - Tq and converting 

millivolts to degrees Fahrenheit. 

Values of the recovery factors given in table I are plotted against 
Reynolds number in figure 3. Included are values of recovery factor ob- 
tained in reference 2 and some previously unpublished data obtained at 
the MCA Lewis laboratory in conjunction with the heat-transfer investi- 
gation presented in reference 4. 

In the turbulent flow region (Reynolds numbers above 3000) the tem- 
perature recovery factor is nearly independent of Reynolds number. The 
average value is 0.88, which agrees with values for flow parallel to 
flat plates . 

For Reynolds numbers below 3000 in the transition region the re- 
covery factor is abruptly reduced to values less than 0.8. The minimum 
and maximum values of recovery factor showed the greatest deviation for 
a Reynolds number of 2900. During this run the static-pressure drop 
along the test section increased slightly, the flow rate was maintained 
constant, and there was an attendant increase in the maximum and minimum 
values of the recovery factor, which indicated the possibility that the 
flow conditions changed from a laminar to a turbulent velocity profile 
near the exit of the test section. ... 

In the laminar flow region the recovery factors increased with de- 
creasing Reynolds numbers. Several check runs were mab.e in this region, 
and the results were not as reproducible as those obtained in the tur- 
bulent flow region. 
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In figure 4 the average friction coefficients "based on the static- 
pressure drop along the test section are plotted against Reynolds num- 
ber. The results indicate that the transition from laminar to turbulent 
flow occurred principally at Reynolds numbers between 2000 and 3000, 
although some deviation from the predicted curve for fully turbulent flow 
occurred at Reynolds numbers up to 6000. 


Analytical Results 

The equations required for calculating recovery factors and heat- 
transfer coefficients were obtained in the ANALYSIS section. 

Effect of wall temperature gradient on recovery factor . - The pres- 
ence of Mach number in equation (l5 ) is due to the fact that the wall 
temperature gradient is not zero. The wall temperature gradient was 
written in terms of Mach number (see eqs. (10), (ll), and (13)) by using 
one-dimensional flow theory, inasmuch as Mach number is a more familiar 
quantity and more readily interpreted. Equations (ll) and (13) or equa- 
tion (B6) indicates that the wall temperature gradient becomes indefi- 
nitely large as the Mach number approaches 1. 

Turbulent recovery factor as a function of Mach number for several 
Reynolds numbers, a Prandtl number of 0.73, and a ratio of eddy ctiffu- 
sivities for heat transfer to momentum transfer of 1 are plotted in fig- 
ure 5. These curves were calculated from equations (15), (16), (17), 
(21), (22), (26), and (27) and the u + against y + plot in figure 1. 
Except for the case M = 0, equation (15) must be solved by Iteration for 
a given value of Tj , inasmuch as T + ’ occurs on both sides of the 

equation. The correct value of T^' was obtained by trial and error. 
That is, it was necessary to try several values of T^ on the right 

side of equation (15) in order to find a value which would agree with 
that given by equation (16) . The curves indicate that the effect of 
wall temperature gradient is negligible at Mach numbers below 0.9 for 
Re = 20,000 and below 0.98 for Re = 390,000. The curve for Hie higher 
Reynolds number is cut off at the point shown because computational dif- 
ficulties were encountered at higher Mach numbers. 

Plotted for comparison are laminar recovery factors obtained by let- 
ting the eddy dlffusivity in the equation for turbulent recovery factor 
equal zero and using the Incompressible parabolic velocity profile for 
laminar flow. These results for laminar flow are probably only qualita- 
tive inasmuch as the one -dimensional approximations made in the analysis 
are not accurate for laminar flow. It is of interest that the effect of 
wall temperature gradient becomes important at much lower Mach numbers 
than for turbulent flow. For M = 0 equation (15) can be integrated, 
and the expression for recovery factor becomes 
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Equation (39) is exact because the one-dimensional approximations are 
not . involved when M = 0. 

Reference 5 states that cp = 2Pr rather than 2Pr - 1. The dif- 
ference is apparently due to the fact that the recovery factor is defined 
in terms of static temperatures (eq. (24)) rather than in terms of total 
temperatures (eq. (23)) in reference 5. As mentioned previously, the two 
definitions would be equivalent only for a = 1 (uniform velocity profile) . 

It should be mentioned that the results for the range of Mach num- 
bers where the recovery factors deviate appreciably from their incompres- 
sible values may be open to question even for the turbulent case, because 
the recovery factor was assumed to be independent of distance along the 
tube. At high Mach numbers the Mach number varies along the tube, so 
that if the recovery factor also varies appreciably with Mach number, it 
will vary along the tube. The analysis may therefore be chiefly useful 
for indicating the range of Mach numbers for which the effect of wall 
temperature gradient is negligible. Fortunately, the results Indicate a 
considerable range of subsonic Mach numbers over which this is true. In 
the remainder of the calculations the effect of wall temperature gradient 
or of Mach number on the recovery factor is neglected. 

Ratio of eddy diffusivlties for heat transfer to momentum transfer 
from recovery factors . - The ratio of eddy diffusivlties for heat trans- 
fer to momentum transfer O' occurs in equation (15) . The ANALYSIS, to- 
gether with experimental data on recovery factors, may therefore provide 
a method for estimating the ratio of the eddy diffusivlties. 

In figure 6 the predicted recovery factors are plotted against Reyn- 
olds number for various values of The curves indicate that <Xj has 

a considerable effect on the recovery factor, the latter decreasing as 
Oj increases. This can be explained by the fact that an Increase In Gl, 
or in the eddy diffusivity for heat transfer, tends to make the static 
temperature more uniform and thus the total temperature less uniform, 
since the two temperatures vary in opposite directions (the static tem- 
perature decreases with distance from the_ wall, whereas the total tem- 
perature increases for recovery factors less than l) . 

In figure 7 the predicted recovery factor curve for an cu of 1.07 
Is compared with the measured values of recovery factors shown in figure 
3. The predicted curve is in reasonable agreement with the data but in- 
dicates a greater variation with Reynolds number than the measured values 
of recovery factor. 

Comparison of the curves for various values of 0^ in figure 6 with 
the measured values indicates that the best agreement would be obtained 
by using a value of that varies from about 1 at a Reynolds number 
of 5000 to 1.09 at a Reynolds number of 400,000. This is the trend 
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that is predicted in reference 8 from a consideration of the effect of 
conduction to an eddy as it moves transversely. It is predicted there 
that ay increases with Peclet number (Pe = RePr) . 

It should be mentioned that the values of O' obtained from recovery 
factors are mean effective values. The local values of <%y may vary with 
radius. More is said about tills in the next section. 

Eddy diffusivity ratio from heat-transfer measurements . - The ratio 
of eddy diffusivities can also be obtained from measurements of heat- 
transfer coefficients and from the ANALYSIS, since equation (35) con- 
tains Gs. Nusselt number as a function of Reynolds number can be calcu- 
lated from equations (35), (36), and (37), or (38), (17), and (25), and 
figures 1 and 8. Hie variation of qj q^ with y* and Tq shown in fig- 
ure 8 was obtained in reference 4. For obtaining Nusselt number as a 
function of Reynolds number, the parameter rj is assigned arbitrary 
values . 

Calculated Nusselt numbers for values of ^ of 1 and 1.07 are 
plotted in figure 9. Included for comparison are experimental data from 
reference 4. The data are in slightly better agreement with the curve 
for an ^ of 1, although the difference between the two curves may be 
within the experimental error of the data. The curves and the data in- 
dicate that, as in the case of recovery factors, the values of Oy at 
the low Reynolds numbers are lower than at the high Reynolds numbers. 

The slightly lower effective values of O' for heat transfer than 
for recovery factors might be explained by assuming that ay increases 
with distance from the wall. In the case of heat transfer the temper- 
ature gradient at the wall is very large compared with gradients at 
other positions, and hence the . effective value of cu should be close 
to the value at the wall. 

In the case of recovery factors, since the temperature gradient Is 
zero at the wall, the effective value of O' should correspond to the 
local value at. some point away from the wall; that is, the effective 
value may be higher. The assumption that Gy increases with distance 
from the wall Is supported by experiments reported in reference 9. 


SUMMARY OF RESULTS 

Measured values of temperature recovery factors for air flowing in 
a tube indicated that: 

1. In the turbulent flow region for Reynolds numbers greater than 
3000, the recovery factor is approximately independent of Reynolds number 
and has an average value of 0.88. 
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2. In the transition region for Reynolds numbers between 2000 and 
3000 the recovery factor is reduced abruptly to a value lower than that 
obtained for the turbulent flow region. 

The following results were obtained from the analysis of recovery 
factors for fully developed flow in a tube: 

1. Wall temperature gradient along the tube has a negligible effect 
on the turbulent recovery factor for Mach numbers less than about 0.9. 

The effect of wall temperature gradient on laminar recovery factors be- 
comes important at much lower Mach numbers . 

2. Comparison of the predicted recovery factors with experimental 
data indicated that the effective value of ratio of eddy diffusivities 
for heat transfer to moment-urn transfer varies from about 1 at a Reyn- 
olds number of 5000 to about 1.09 at a Reynolds number of 400,000. 

The effective values of eddy diffusivity ratio appeared to be slightly 
less for heat-transfer coefficients. 

3. The analysis substantiated the experimental result that correla- 
tions for heat transfer without frictional heating can be applied to 
frictional heating by basing the heat-transfer coefficient on the differ- 
ence between the wall temperature and the adiabatic wall temperature 
(for the same bulk temperature) , rather than on the difference between 
the wall and the fluid bulk temperatures . 


Lewis Flight Propulsion Laboratory 

Rational Advisory Committee for Aeronautics 
Cleveland, Ohio, July 22, 1958 
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APPENDIX A 


DERIVATION OF TURBULENT ENERGY EQUATION 

The method of derivation given here is essentially the same as that 
used in reference 10 for flow over a flat plate. The equations of mo- 
mentum, continuity, and energy used herein for symmetrical flow in a tube 
can be written, respectively, as 


pu 


3u , chi 1 b ( 3u\ 

33E + pu r FF = x 3? " 


3(pu) . 1 

X! H 


"S 


3F 


= o 


(Al) 

(A2) 


puc T 


H + t? - h ^ ( rk I?) + r (If) + u li (A3} 


Time derivatives are neglected in these equations, as flow is assumed 
statistically steady. Multiplying equation (Al) through by u and add- 
ing the result to equation (A3) give 


pu 


b( u 2 \ 3 / u 2 \ 1 3 f . 3t 3u\ /. v 

S (V + xj + pu r S (V + T) = ;5|' k S tI,B 5) < A4) 


where Cp is considered constant. If the instantaneous quantities in 
equation (A4) are replaced by the sum of their time average and fluctuat- 
ing components, and the equation of continuity is used, the following is 
obtained : 

?5 s (°p 5 + V) + w*r + p 7 ^:) ^ (y + 1 $ f r ( 5 1? + 


(y + v) + + p, y 3? (v *l)'75 [ r ( k 

H - P= p “i* 1 - p“ “’“ijl 


(A5) 


where the boundary-layer assumptions have been used and the fluctuations 
of the physical properties are assumed small compared with their mean 
values. But because 
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equation (A5) becomes 

( c p t + v) + ^ s (y + v) = § -I [ r ( 5 if + >“ if - 

pe p u’t> - pu u’u^yj (A6) 

For fully developed flow there is no mean radial mass flow, or, pu^. is 

zero. If, in addition, r is replaced by (r Q - y), u? *>y ~ v > and tlie 

left side of equation (A6) is written in terms of total temperature, it 
becomes - _ ■ y, 

(r 0 - y)c p pu |jj = ^ [(r 0 - y) ^ - P=p ~ ^ - P ; ^jj < A7 > 
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APPENDIX B 


CALCULATION OP dT 0 /dx 


The tenqperature recovery factor is given in terms of total hulk 
temperature as 


cp = 1 - 



(Bl) 


For fully developed flow the recovery factor is independent of x. Dif- 
ferentiating equation (Bl) and using equation (7) give 


<™0 = fob - g p) du b 

dx u^/ 2 dx 


(B2) 


In the present development one-dimensional flow is assumed; that is, the 
fact that the velocity and temperature vary with radial position is ne- 
glected. This is a reasonable assumption for turbulent flow. By using 
the continuity relation w = P^u^A, equation (B2) can be written as 


2(T b - T q ) 
dx ~ dx 

In order to obtain dp^/dx in terms of known quantities, the total- 

pressure drop is written as the sum of the momentum and the friction 
pressure drops, or 


dp 2 ^b 2f 2 

dx ~ U B dx " D 


(B4) 


Differentiating the perfect gas law p = p^Rt^ and the definition 
of total bulk temperature for one -dimensional flow (a = l) 


Ux 


T b S *b + 25- 


. 2c^ 

a* p 

and using equation (7) and continuity give on substitution into (B4) 

2 , 


fPb 

dx 


ffi b u b/ r O 

4 - Ru I/ c p - Rt b 


(B5) 
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Substituting equation (B5) into (B3) and using the perfect gas relation 
hp = Ry/ (y - l) and the definition of Mach number give 


dl 0 2(T b -T 0 )f yM 2 


ax 


(B6) 


L 0 1 - M 

or, using the definitions of T^ , and f, the resulting equation is 

,2 


ZCf ' wi 

M 


1 - 1(2 


(B7) 
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temperature 
of fluid 
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Figure 4. - Correlation of isothermal average friction factor with Reynolds number evaluated at bull total 
temperature of fluid. 
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Figure 5. - Variation; of recovery factor vltb Mach nuriber at Erandtl nunlber of 0.75 and several Reynolds munbers. 
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Figure 8. - Variation of beat transfer per unit area across tube. 
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